Localization and Pattern Formation 
in BBGKY Hierarchy 



Antonina N. Fedorova, Michael G. Zeitlin 



Mathematical Methods in Mechanics Group 
Institute of Problems of Mechanical Engineering (IPME RAS) 
Russian Academy of Sciences 
Russia, 199178, St. Petersburg, V.O., Bolshoj pr., 61 
zeitlin@math.ipme.ru, anton@math.ipme.ru 
http://www.ipme.ru/zeitlin.html 
http://www.ipme.nw.ru/zeitlin.html 



Abstract 



A fast and efficient numerical-analytical approach is proposed for mod- 
eling complex behaviour in the BBGKY hierarchy of kinetic equations. 
Numerical modeling shows the creation of various internal structures 
from localized modes, which are related to the localized or chaotic type 
of behaviour and the corresponding patterns (waveletons) formation. 
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We consider the application of a new numerical/analytical technique based on local nonlinear harmonic analysis approach 
(multiscale decomposition) for the description of complex (non-equilibrium) behaviour of the classical statistical ensem- 
bles, considered in the framework of the general BBGKY hierarchy (kinetics equations). We restrict ourselves to the ratio- 
nal/polynomial type of nonlinearities (with respect to the set of all dynamical variables) that allows to use our results from 
[1],[2], which are based on the so called multiresolution framework [3] and the variational formulation of initial nonlinear 
(pseudodifferential) problems. Our main goals are an attempt of classification and construction of a possible zoo of nontrivial 
(meta) stable states: high-localized (nonlinear) eigenmodes, complex (chaotic-like or entangled) patterns, localized (stable) 
patterns (waveletons). The last case is a good image for plasma modeling in fusion state (energy confinement). It should be 
noted that the class of smoothness (related at least to the appearance of chaotic/fractal-like type of behaviour) of the proper 
functional space under consideration plays a key role in the following [2]. Let M be the phase space of an ensemble of N 
particles (dimM = 6N) with coordinates x t = (qi,Pi), i = l,—,N, q { = (qj,qf,qf) G R 3 , p { = (pj,p 2 ,p 3 ) G 
R 3 i q = (<Zi> • • • j ?Jv) £ R 3N ■ Individual and collective measures are: pi = dx-i = clqi&pi, p — Yli=i Mi- O ur construc- 
tions can be applied to the following general Hamiltonians: Hn = J2iLi (l^l + ^(<?)) +J2i<i<j<N Uij(q%,qj), where the 
potentials Ui(q) = Ui(q\, . . . , qjy) and Uij(qi, qj) are restricted to rational functions of the coordinates. Let L s and Lij be 
the Liouvillean operators (vector fields) and Fn(x\, . . . , xn; t) be the hierarchy of TV-particle distribution function, satisfying 
the standard BBGKY-hierarchy (V is the volume): 

OF 1 f x— -\ 

' • L S F S = — / dfjb a+ i ^2 L iiS+ xF s+ i (1) 

i— 1 



dt V s 



Our key point in the following consideration is the proper nonperturbative generalization of the perturbative multiscale ap- 
proach of Bogolyubov. The infinite hierarchy of distribution functions satisfying system (1) in the thermodynamical limit 
is: F = {F , Fx(xx;t), . . . , F N (xx, . . . ,x N ;t), . . .}, where F p (x 1 ,...,x p ;t) G H p , H° = R, H p = L 2 (R 6p ) (or any 
different proper functional space), F G H°° — H° © H 1 © ■ • • © H p © . . . with the natural Fock space like norm (guaranteeing 
the positivity of the full measure): 

(F,F) =F? + J2 I F?(x u ...,x i ;t)f[ t i i . (2) 



First of all we consider F = F(t) as a function of time only, F G L (R), via multiresolution decomposition which nat- 
urally and efficiently introduces the infinite sequence of the underlying hidden scales. Because the affine group of trans- 
lations and dilations generates multiresolution approach, this method resembles the action of a microscope. We consider 
a multiresolution decomposition of L 2 (R) [3] (of course, we may consider any different and proper for some particu- 
lar case functional space) which is a sequence of increasing closed subspaces Vj G L 2 (R) (subspaces for modes with 
fixed dilation value): ...V-2 C V-i C Vo C V\ C V% C .... The closed subspace Vj(j G Z) corresponds to the 
level j of resolution, or to the scale j and satisfies the following properties: let Wj be the orthonormal complement of 
Vj with respect to Vj+\\ Vj+\ = VjQWj. Then we have the following decomposition: {F(t)} = ©_ 00< j <00 Wj = 

oo 

V c Wj in case when V c is the coarsest scale of resolution. The subgroup of translations generates a basis for the 

3=0 

fixed scale number: apan keZ {2 J ^ 2l i(2 : >t — k)} — Wj. The whole basis is generated by action of the full affine group: 
sparing z j eZ {2 J / 2 ~$>(2H — k)} — span fe j^z{^i,k) — {F{t)}. In multidimensional case we may consider polynomial 
tensor bases, e.g. in n — 2 case we may use the rectangle lattice of scales and one-dimensional wavelet decomposition: 
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f(xi, X2) — Y^i e-j fe(/> ^i,t ® ^j,k)^j,e ® ^j.fcOci, X2), where the basis functions ® depend on two scales 2 4 
and 2 _J . We obtain our multiscale/multiresolution representations (formulae (3) below) via the variational wavelet approach 
for the following formal representation of the BBGKY system (1) (or its finite-dimensional nonlinear approximation for the 
rt-particle distribution functions) with the corresponding obvious constraints on the distribution functions. Let L be an arbi- 
trary (non)linear differential/integral operator with matrix dimension d (finite or infinite), which acts on some set of functions 

from L 2 (S!®"): * = \P(t, x\, x 2 , ■ ■ ■ ) = ($! 1 (t,Xi,X2, ...),..., ^ d (t, Xi, x 2 , ... )J, Xi 6 SI C R 6 , n is the number of 
particles: L* = L(Q,t,Xi)^(t,Xi) = 0, Q = Qd ,d u d 2 ,...(t, Xi, x 2 , d/dt,d/dxi, d/dx 2 , ...,Jnk) = ^tof*/^ ■'■'■=! 
Qioiii 2 — (*> x \i x 2, ■ ■ ■ ) (w) (air) ( sfj) ' • • / t ik ~ ^ et us cons idei" the N mode approximation for the solution as the 
following ansatz: ^ N {t,xi,x 2 , ■ ■ ■) = Y,i ,i u i 2 ,-=i a i iyi2..Ai () ®C i2 . .. (t, x\, x 2 , . . . ). We shall determine the ex- 

pansion coefficients from the following conditions: £^ q klk2 = J (m> N )Ak (t)Bk 1 (xi)Ck 2 (x 2 )dtdx±dx 2 ■ • ■ = 0. Thus, 
we have exactly dN n algebraical equations for dN n unknowns OtQ,^,.... This variational approach reduces the initial problem 
to the problem of solution of functional equations at the first stage and some algebraical problems at the second. The solution 
is parametrized by the solutions of two sets of reduced algebraical problems, one is linear or nonlinear (depending on the 
structure of the operator L) and the rest are linear problems related to the computation of the coefficients of the algebraic 
equations. The solution of the equations (1) has the following multiscale decomposition via high-localized eigenmodes 

F(t,xi,x 2 , ■ ■ ■) = a l3 U l ®V'{t,x ll ...) 1 (3) 

y ]{t) = V ^o W{t) + £ V j {uit)j ^ ^ ^ w[xs) = V ^o W{xs) + J2 UUk^Xs), k* m ~ 2™, 
1>N m>M 

which corresponds to the full multiresolution expansion in all underlying time/space scales. The formulae (3) give the expan- 
sion into a slow and fast oscillating parts. So, we may move from the coarse scales of resolution to the finest ones for obtaining 
more detailed information about the dynamical process. In this way one obtains contributions to the full solution from each 
scale of resolution or each time/space scale or from each nonlinear eigenmode. It should be noted that such representations 
give the best possible localization properties in the corresponding (phase)space/time coordinates. Formulae (3) do not use 
perturbation techniques or linearization procedures. Numerical calculations are based on compactly supported wavelets and 
related wavelet families [3] and on evaluation of the accuracy on the level N of the corresponding cut-off of the full system (1) 
regarding norm (2): ||F Ar+1 — F \\ < e. Fig. 1 demonstrates the appearance of localized (meta) stable pattern (waveleton), 
which can be considered, e.g., as a model for fusion state in plasma. 
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Fig. 1 Localized waveleton pattern. 
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